Viscous to Inertial Crossover in Liquid Drop Coalescence 
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Using an electrical method and high-speed imaging we probe drop coalescence down to 10 ns after 
the drops touch. By varying the liquid viscosity over two decades, we conclude that at sufficiently 
low approach velocity where deformation is not present, the drops coalesce with an unexpectedly 
late crossover time between a regime dominated by viscous and one dominated by inertial effects. 
We argue that the late crossover, not accounted for in the theory, can be explained by an appropriate 
choice of length-scales present in the flow geometry. 
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Typically, it is a simple exercise to estimate when a 
fluid flow will switch between a viscous and an inertially 
dominated regime. The dimensionless Reynolds number 
serves this purpose: by identifying the characteristic ve- 
locity (U) and length (L) one expects crossover behavior 
when Re = pUL/p w 1 where p and p are the liquid 
density and viscosity respectively. In this way, Eggers et 
al. [EHI compute a Reynolds number for liquid drop co- 
alescence, which predicts that for the case of salt water, 
viscous forces give way to inertial ones just 0.7 ns after 
the drops touch. 

High-speed imaging experiments |3|-[6j have observed 
the coalescence of water drops at speeds up to 10 6 frames 
per second, but limited spatial and temporal resolution 
prevented these studies from observing the initial regime 
dominated by viscous effects. In this paper, we use an 
electrical method 0, H[ to observe salt water coalescence 
down to 10 ns after the drops touch. Our measure- 
ments are conducted at slow enough approach velocities 
where no drop deformation occurs. Followingthe widely 
accepted Reynolds number for coalescence [3-0? @L we 
would expect to see only the inertial regime. However, 
we observe viscous effects until 2 ps after contact which 
is more than 3 decades longer than predicted. 

We argue that the source of this discrepancy is that the 
previously used Reynolds number is based on an incor- 
rect length-scale. Our detailed measurements covering 
two decades in liquid viscosity suggest a new picture of 
the crossover; the correct Reynolds number for drop coa- 
lescence is based on a length-scale, not fully appreciated 
previously, that reflects the dominant flows. 

Experiment — As illustrated in Fig. Ufa), liquid drops 
are formed on two vertically aligned teflon nozzles of ra- 
dius A = 2 mm. One drop is fixed while the other drop 
is slowly grown with a syringe pump until they coalesce. 
Our experiments were performed at ambient air pressure. 
The intervening gas layer between two colliding drops can 
distort the drops and delay their coalescence [l0[ . Previ- 
ous experiments @, Q suggest that such distortion may 
be present for approach velocities, U app , as low as 10 -4 
m/s. To this end, we have measured the effect of the 
ambient gas. The dynamics reported in this paper are in 



a low U app regime (U app <8.0xl0 _5 m/s) where the gas 
does not disturb the initiation of coalescence. 

We follow the electrical method developed by Case et 
al @, H| to isolate the time- varying complex impedance, 
ZcRi of two liquid hemispheres as they coalesce. A high- 
frequency (0.6 < / < 10 MHz) AC voltage, V$ n , is ap- 
plied across the drops at low amplitude (V* n < 1 V). By 
measuring two voltages, V\ and V2 shown in Fig. Ufa), 
we extract Zcr, which we model as a time- varying resis- 
tor (Rcr) and capacitor (Ccr) in parallel. By applying 
an additional DC offset voltage, we can determine that 
the electric fields do not affect the measurement of the 
coalescence. A sharp feature in the phase of V2 at the in- 
stant the drops touch allows us to determine the moment 
of contact, to, to within 1//. 

We calculate the conversion between Rcr and the neck 
radius r using the electrostatics calculation package ES- 
tat (FieldCo) p], [8[ • We compared the calculation of two 
possible bridge geometries while fixing the electrical po- 
tential at the nozzles. The results agree with each other. 
This implies that the minimum bridge radius, r, is the 
single geometrical feature controlling the resistance. We 
find an excellent fit to: 



Rcr = 2/^ar + Ro 



(1) 



from r = out to r = A/3, where a is the electrical con- 
ductivity of the fluid, £ = 3.62±0.05 is a fitting parameter 
obtained from the simulations, and Ro = I /an A. 

According to this conversion, the measured quantity of 
interest is Rcr — Ro, which we show for aqueous NaCl 
coalescence in Fig.[TJb) as a function of r = t — to- In the 
inset, we convert this measurement to a bridge radius, r, 
and show that it agrees with high-speed imaging results. 
This comparison demonstrates not only the quantitative 
accuracy of the electric method but also its superior dy- 
namic range as compared to the optical techniques. 

Comparison to theory — Coalescence begins in a viscous 
regime, where surface tension, 7, is balanced by viscous 
forces. Based on the analytic solution in two dimensions 
[Tlj . it is predicted that for three-dimensional drops [l[: 
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FIG. 1: (a) Coalescence cell and measurement circuit. We ap- 
ply an AC voltage, Vi n (Hewlett-Packard, HP3325A), across 
known circuit elements (Rk and Ck) and coalescing liquid 
drops. We read voltages V\ and V2 into Labview (National 
Instruments) to calculate the cell impedance (dotted box). 
Zcr is from the coalescence region (dashed box); Z t and Zb 
are due to the top and bottom nozzles; C p is a small cell ca- 
pacitance in parallel with Zcr] and Ri n = 50f2 is the input 
impedance of a high-speed digitizer (National Instruments, 
NI PCI-5105). The nozzles are brought together to measure 
Z t + Zb. (b) Mean value of Rcr — Ro for six coalescences 
of aqueous NaCl (7 88 mN/m, p = 1180 kg/m 3 , /1 = 1.9 
mPa s). Vertical error bars are the standard deviation of the 
points in each logarithmically spaced bin. Horizontal error 
bars are ±1/(2/). Asymptotic behavior is consistent with 
1.3 x 10" 3 r _1 at early times (dotted line), and 0.90r" 1/2 at 
late times (dashed line). Inset: Measurements of the bridge 
radius from the electric method (o) probe much earlier than 
high-speed imaging (x), but extend just as far to long times. 



tension, a scaling argument [1] produces: 



' inviscic 
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(4) 



Numerical simulations reproduce this scaling 0, [H, 13 ] 
and give Do ^1.62 [2|. High-speed imaging experiments 
(^l-lil Il2| also observe this scaling regime. 

For 2 mm drops of aqueous NaCl solution in air, we 
therefore expect Rcr — Ro = 4.5 x 10 _4 r _1 at early 
times, and Rcr — Ro = l.lr -1 / 2 at late times (using 
eqns. [TJEJ andHJ. As shown in Fig. H^b), our measured 
Rcr — Ro is in agreement with both of these predicted 
asymptotic scalings. 

However there is a glaring discrepancy. Theory pre- 
dicts that the crossover time between these regimes, r c , 
is roughly 0.7 ns, whereas we see r c « 2 /is, which is more 
than 3 decades later than expected. We investigate this 
discrepancy by varying the liquid viscosity. 

Varying the liquid viscosity — Our liquids, mixtures of 
glycerol and water, were saturated with NaCl to make 
them electrically conductive. We measured the viscos- 
ity, surface tension, density, and electrical conductivity 
of each fluid. By varying the glycerol content, the vis- 
cosity could be varied over two decades (from 1.9 mPa 
s to 230 mPa s) while the surface tension and density 
remained nearly constant, changing by only a factor of 
1.6 and 1.04, respectively. 

The predicted viscous scaling (eqn. [2]) with the loga- 
rithmic correction (eqn. [3j) is an asymptotic result that 
has no free parameters and is predicted to apply for 
r/A < 0.03 [l|, which is difficult for optical methods to 
probe. Fig. [2fa) shows our data from electrical measure- 
ments that extend 2 decades below 0.03A We compare 
this data with Co given by eqn. [3] (dashed lines) and with 
Co = 1, i.e., linear expansion at the capillary speed, 
(solid lines). The data are better fit by the latter; we find 
no evidence for the predicted logarithmic corrections. 

The inset to Fig. Efb) shows r(r) for ten viscosities, 
ranging from 1.9 mPa s to 230 mPa s. As shown in Fig. 
[21(b), the data collapses cleanly onto itself when the axes 
are rescaled: r — » r/r c and r — >• r/r c where r c and r c 
are free parameters at each viscosity to produce the best 
collapse. The entire data in the master curve can be well 
fit with the simple interpolation: 



where Co is calculated to be: 
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We expect Co to be nearly unity over our measurement 
range. High-speed imaging experiments 0, 0, El that 
corroborate eqn. [2] measure prefactors of order unity, but 
are not sensitive to the logarithmic corrections. 

In the other limit, where inertial forces balance surface 
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This collapse determines the coefficients for the early and 
late-time scaling. We compare these coefficients to the 
predicted values in Fig. [3j our measurements of Co are of 
order 1 across the entire range of viscosity, and Do is in 
good agreement with the predicted value of 1.62 [2|. 

At low viscosity, there is a small departure from unity 
in our measurement of Cq. A possible cause is the non- 
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FIG. 2: (color online) Bridge radius, r, versus r for glycerol- 
water-NaCl mixtures of various viscosities. At each viscosity, 
measurements from 5 or more runs are logarithmically binned 
and averaged, (a) At two viscosities, data for r(r) are com- 
pared with eqn. [2j with logarithmic corrections given by eqn. 
[3] (dashed lines) and with Co = 1 (solid lines). The data fit 
better to Co = 1. (b) Inset: r(r) for viscosities ranging from 
1.9 to 230 mPa s. Main: Data are collapsed by rescaling 
the horizontal and vertical axes by r c and r c for each viscos- 
ity. Limiting behavior is proportional to r/r c (dotted line) 
at early times, and y/r /r c (dashed line) at late times. The 
collapsed data are well described by eqn. [5] (solid line). 

vanishing dynamic viscosity ratio between the surround- 
ing air and our fluids. This would be larger at low liquid 
viscosity, consistent with the data. These effects should 
be negligible for p > 82 mPa s, the lower viscosity curve 
in Fig.|2ja), where the ratio is 4 x 10 4 . We also note that 
if the logarithmic corrections of eqn. [3] hold, there should 
be small deviations from the master curve, Fig. [2fb). 
These corrections are difficult to access experimentally 
but would be more pronounced as one goes to smaller r. 

The crossover time, r c , as a function of viscosity, /i, is 
shown in Fig. [3jc). The data are fit well by a quadratic 
dependence on p (solid line). Clearly, the accepted for- 
mula for r c oc p 3 (dashed line) does not agree with the 
data. This suggests that the conventional Reynolds num- 
ber for coalescence, Re = pj 2 r/p 3 ^ is wrong. This esti- 
mate is based on using r and dr/dr as the dominant 
length- and velocity- scales. We argue that these are not 
the right choices. Instead, we suggest that the dominant 
flows occur over a much smaller length: the neck height. 

Reynolds number for coalescence — We can equally well 
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FIG. 3: (a and b) Measured dimensionless scaling- law pref- 
actors, Co and Do, versus viscosity. In (a), the dashed line 
is Co = 1. In (b), the dashed line is Do — 1.62, the value 
obtained from simulation [2j]. Error bars show the range in 
which r c and r c can be varied without affecting the quality 
of the collapse, (c) Viscous to inertial crossover times versus 
viscosity, obtained via the collapse in Fig. [2] The solid line, 
t c = 0.3/i 2 , is based on our proposed Reynolds number (i.e. 
eqn. [6] with p 1200 kg/m 3 , 7 = 65 mN/m, D 1.62) and 
fits the data well. The dashed line shows r c = p 3 / pj 2 = 0.2p 3 
obtained from the Reynolds number proposed in the literature 
0,0. Clearly this is a poor fit to the data. 

prescribe a Reynolds number coming from either the vis- 
cous or the inertial regime. In the inertial regime, the 
neck height (the spacing between the drops just outside 
the fluid neck) is r 2 /A [2]. As illustrated in Fig.SJ liquid 
from each drop moves to fill half of the gap, L « r 2 /2A. 
We assume that the flow occurs in an annular region hav- 
ing a width equal to the gap size. Due to volume consid- 
erations (see Fig. H]) , the flow speed is proportional to the 
interfacial velocity, U « \drjdr. Using eqn. [4] for r, we 
get L « iDlin/pAf^r and U * \D^A/ pY^r' 1 ' 2 . 

Using these characteristic scales, the Reynolds number 
is: Re = pUL/p w £>o(/?7 3 jAf^r 1 ' 2 /8p. The crossover 
time, r c , is obtained by setting Re = 1: 

For our fluids, we find r c « 0.3/i 2 , which is in good agree- 
ment with our data (Fig. EJc)). Hence, our proposed 
Reynolds number gives crossover times consistent with 
our experiments ranging over two decades in viscosity 
[tH - The only place where the two predictions (r c ex p 2 
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FIG. 4: Length- and velocity- scales for coalescence, (a) In 
the inertial regime, a fluid bridge of radius r has height r 2 /A. 
(b) As the neck expands radially, fluid fills in the neighboring 
region with a flow, which occurs on a scale set by half the 
gap size, L w r 2 /2A, that is essentially axial in direction 
(vertical arrows). The speed of this flow, U, can be related to 
dr/dr. The volume sandwiched between the drops out to a 
radius r (shaded region) is given by V = ^r 4 . As the bridge 
expands, dV/dr = ^-r 3 dr /dr. This fluid is supplied through 
two annular regions of radius r and width r 2 /A, above and 
below the gap (dotted lines), of total area -^f-r 3 . Setting the 
volume flow rate, ^r 3 [7, equal to dV/dr, we get U « \ dr/dr. 

versus r c ex fi 3 ) give similar results is at high liquid vis- 
cosity, which is where two previous crossover measure- 
ments had been reported 0, [H[ . 

This argument can be recast in terms of a crossover 
length. Using eqns. |4]and[6j we find Re « 1 corresponds 
to when r 2 / A « M 2 /P7 = l v , where l v is the viscous 
length-scale. In other words, coalescence proceeds in the 
viscous regime until the gap between the drops becomes 
as large as the viscous length-scale of the fluid. Therefore, 
the relevant length scale to which l v should be compared 
is not the bridge radius, r (as in p], 0, @-[9|, [HI), but 
rather the bridge height, r 2 /A, which is much smaller. 
Finally, two- and three-dimensional coalescence are ex- 
pected to be equivalent to leading order Thus the 
above argument also applies to 2D coalescence and ex- 
plains an unexpectedly large crossover radius measured 
in the coalescence of liquid lenses [15[ . 

Conclusion — We have probed coalescence down to 10 
ns after the drops touch, providing a detailed study of 
the viscous to inertial crossover dynamics for liquid drop 
coalescence. The surprisingly late crossover we observe 
supports a new picture where a length scale drastically 
smaller than the bridge radius controls the flow. In the 
viscous regime, the data are better fit with a constant 
expansion velocity than the form predicted to have log- 
arithmic corrections. One remarkable outcome is that 
our data, over a wide range of viscosity and time, can be 



rescaled onto the master curve in Fig. O This includes 
a very long crossover region between the viscous and in- 
ertial regimes, where there is no quantitative theoretical 
work describing the bridge dynamics. 
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